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We have re-examined the Andreev-Lifshitz theory of supersohds. This theory implicitly neglects 
uniform bulk processes that change the vacancy number, and assumes an internal pressure P in 
addition to lattice stress Xit- Each of P and Xik takes up a part of an external, or applied, pressure 
Pa (necessary for solid *He). The theory gives four pairs of propagating elastic modes, of which 
one corresponds to a fourth-sound mode, and a single diffusive mode, which has not been analyzed 
previously. The diffusive mode has three distinct velocities, with the superfluid velocity much 
larger than the normal fluid velocity, which in turn is much larger than the lattice velocity. The 
mode structure depends on the relative values of certain kinetic coefBcients and thermodynamic 
derivatives. We consider pressurization experiments in solid ^He at low temperatures in light of this 
diffusion mode and a previous analysis of modes in a normal solid with no superfluid component. 

PACS numbers: 67.80. bd, 05.70.Ln 



I. INTRODUCTION 

Since the late 1960's there have been theoretical sug- 
gestions that solids might display flow behavior similar to 
what is found in superfluids>i"— For that reason there has 
been a great deal of interest in solid ^He as a candidate to 
be a supersolid.^ The first experimental indication of su- 
perflow was the appearance of a non-classical moment of 
inertia (NCRI), first observed by Chan's group, since con- 
firmed by many other laboratories, and strongly linked 
to disorder!^"— In addition, the shear modulus shows 
anomalous behaviorji^ although not enough to explain 
the NCRI experimentsjii Non-NCRI superflow has been 
searched for but not observed fiS. Evidence is growing that 
restricts the possible temperature range over which su- 
persolidity can occur. ^- Moreover, for NCRI experiments 
with rim velocity v at temperature T, the observed hys- 
teresis in V — T space suggests multiple apparent phase 
transitions 

A recent experiment on a pancake-shaped sample, 
where a pressure change is applied to one side, finds 
an exponential decay with time of the pressure response 
on the opposite side^^S, The response is slower at lower 
temperatures, rather than saturating as for a quantum 
transition, perhaps an indication that the system is not 
supersolid. We have recently studied the lattice diffusion 
mode of a normal solid (see Refl23|). based on equations 
obtained by eliminating the superfluid velocity Vg from 
the theory of Andreev and Lifshitz.^ We obtained both 
the diffusion constant and the eigenmode structure (by 
which we mean the ratios of the deviations from equi- 
librium of various thermodynamic quantities) for a solid 
under an externally applied pressure Pa (necessary to so- 
lidify ^He, even at T = 0). 

Whereas an ordinary solid has eight degrees of 
freedom)^ the addition of Vg (the gradient of a phase) 
gives a supersolid nine degrees of freedom. For a plane 
wave, where ki is the wavevector with magnitude k) the 
degrees of freedom are given by two scalar thermody- 
namic quantities (which can be taken to be the mass den- 



sity p and the entropy density s), the lattice vector Ui, the 
normal fluid velocity vector Vni, and Vs — kiVsi/k, where 
Vsi is the superfluid velocity. The scalar quantity Vg has 
been defined because Vsi is expressible as the gradient of 
a phase 0. The total momentum density is thus given by 
ffj = PnVni + PsVsi, where p„ and ps are the respective 
densities of the normal and superfluid components. (In 
principle, both p„ and ps are tensors, but calculations 
for hep *He indicate that they are nearly isotropic— i^) 
The nine degrees of freedom imply there are nine nor- 
mal modes. For a uniform inflnite system these modes 
are: four pairs of propagating elastic waves (previously 
studied for both zeroi^^^^— and nonzero^^ Pa), with fre- 
quency uj ^ k; and a diffusive mode, with w ~ ik'^ (whose 
structure has not been previously studied). 

The present work studies this diffusive mode that oc- 
curs in a supersolid when Vg is included. We employ a 
variation on the notation of ReflU, which gives a more 
explicit derivation of the equations of motion than does 
Ref,h, and extends Ref, l] to include nonlinear terms 
Ref.l and Ref, 26 implicitly assume that uniform vacancy- 
number-changing bulk processes are negligible. 

The Andreev-Lifshitz theory is remarkable in that it 
assigns an internal pressure P, in addition to lattice stress 
Aifc, to a supersolid, in order to continuously go to the 
superfluid limit. Each of P and Xik take up part of Pa- 
Ref. m finds for a solid, by thermodynamic considera- 
tions, the dependence of P on P^. The consequences of 
distinct P, Xjh, and Pa had not previously been consid- 
ered. Ref. [23 calculates the effect of Pa on the prop- 
agating elastic and diffusive modes of an ordinary solid. 
Ref. |29| calculates the effect of Pa on the propagating elas- 
tic modes of a supersolid, as well as the efficiency with 
which a heater or a transducer generates these modes. 
The present work considers the effect of Pa on the diffu- 
sive mode of a supersolid. 

As for the lattice diffusion mode for the normal solid, 
the diffusive mode for the supersolid is characterized not 
by the diffusion of a single thermodynamic variable, but 
by specific amounts of each, determined by the eigen- 
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mode structure. A dissipative term in the equation of 
motion for the lattice displacement permits the lattice 
velocity to differ from u„. We find the relationships be- 
tween the normal, superfluid, and lattice velocities in this 
mode. From the lattice velocity one can obtain the lattice 
displacement and lattice strain deviation. Because the 
mode is diffusive, the rate of change of the momentum 
density, and thus the total stress deviation, are nearly 
zero, so lattice stress deviations must be canceled by 
an opposing pressure deviation P', thus determining P' . 
Again because the mode is diffusive, the rate of change 
of Vs, and thus the total chemical potential deviation, 
are nearly zero, so the P' must be compensated by an 
opposing contribution due to a temperature deviation 
T'. This diffusion mode is therefore characterized by its 
diffusion constant and specific ratios of the normal and 
superfluid velocities, and the temperature and pressure, 
relative to the lattice velocity. In practice we use the 
entropy and mass densities rather than temperature and 
pressure. The theory permits vacancies to diffuse but 
there are no bulk sources or sinks for them. 

Section HI] gives the AL supersolid theory in our nota- 
tion. Section IIIII derives the normal modes for the su- 
persolid. Section IIVI provides a summary. Appendix |X] 
estimates the sizes of several quantities relevant to the 
diffusive mode. 



II. ANDREEV-LIFSHITZ SUPERSOLID 

In what follows we employ the primary quantities 
energy density e, lattice displacement Ui, and non- 
symmetrized strain Wik = diUk- 

A. Thermodynamics 

The thermodynamic equations for a supersolid are 
given by 

de = Tds + Xikdwik + ^.dp + w„ • dg + js ■ dvg, (f ) 

e = -P + Ts + KkWik + fJ.p + Vn ■ g + js ■ Vs, (2) 

= -dP + sdT + WikdXik + pdfj. + g ■ dvn + Vs ■ djs- 

(3) 

Here A^^ is an elastic tensor density (with the same units 
as pressure P), p is the chemical potential (with units of 
velocity squared), 

g ^ PnVn+ PsVs (4) 

is the momentum density, and 

js = g- pVn = Ps{Vs - V,i)- (5) 

is a momentum density defined so that de ^ g ■ d{Sv) 
under a Galilean boost Sv. Since (g, VmVs) are all vectors 
under Galilean boosts, we deduce that Pn + Ps = P- 



We find it convenient to define 

jn=pVn, (6) 

so that 

g^3n+js- (7) 

Unlike js, the quantity jVi is a momentum density both 
in units and in its properties under Galilean boosts. 

B. Dynamics 

The linearized equations of motion relevant to obtain- 
ing the normal modes, considering only the independent 
variables s, u^, p, and Vsi, are 



5ts + a,/, = (i?>0) (8) 

dtu, = U^, (9) 

dtp + d,g, = 0, (10) 

^tg^ + ^kn,k = Q, (11) 

dtvs, + d^e = 0. (12) 

Here, the fluxes /i, liik, gt, 0, and the "source" Ui are 
given by 

f^ = sv^,-^d,T~^diX^„ (13) 

U^ = + ^djT + p.jdiXij , (14) 

Ilifc = [PSik - Xki) ~ ViklmdrnVnl " Cikdljsl, (15) 

S = P - CikdkVni ~ Xdkjsk, (16) 

PVm+jst, (17) 



and we take i? « 0, as it is second order in deviations. 
AL use both aik ~ —^ik (a notation we employ below) 
and ji = gi. The term in ()14p proportional to /3y al- 
lows the lattice velocity iii to differ from the velocity u„j 
associated with mass flow. 

Recall that a diffusion constant D is proportional to 
a characteristic velocity times a characteristic mean-free 
path, so it has units of m^/sec. In terms of a D, the 
dissipative coefficients have the following units: Hij has 
units of s times D; aij has units of D; l3ij has units of 
inverse pressure times D; rjikim has units of p times D] 
Qik has units of D; and x has units of inverse density 
times D. 



III. NORMAL MODES IN A SUPERSOLID 

As noted earlier, this system has nine variables: s, 
p, Ui, Vni, and Vg- With deviations from equilibrium 
denoted by primes, we use the nine variables s', p' , u\, 

gl^ PnVn',+ PsVs'i, (18) 
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and 



(19) 



As noted above, there correspondingly are nine normal 
modes. For an infinite system we assume a disturbance 
of the form exp[z(A: • r — ujt)], where the real wavevector 
k is considered to be known, but w is unknown. For the 
disturbance to decay in time, Im{uj) < 0. We find that 
six modes come in three degenerate pairs, with g'l and 
u'^ strongly coupled, and correspond to ordinary elastic- 
ity. Two other modes also form a degenerate pair, cor- 
responding to fourth sound, with the superfluid compo- 
nent in motion and the normal component essentially at 
resti^i^i^ The ninth and final mode is diffusive, with 
and v'g in opposing directions, and nearly constant 
chemical potential and stress. 

We consider the (off-diagonal) temperature-lattice 
transport coefficient a^- = 0, and set to zero the ther- 
mal expansion coefficient. We also neglect the viscosities 
Vikim, Cik, and X, which to lowest order do not contribute 
to the modes."^-S We consider an isotropic solid, for which 
Ktj = k5ij and /3y = pSij. 

Unless otherwise specified, thermodynamic derivatives 
with respect to p, s, or Wik are taken with the other two 
variables held constant. 



A. Elastic Modes 

The elastic modes are obtained by neglecting dissipa- 
tive and nonlinear terms in ([5|)- p^ . Although the elas- 
tic modes of a supersolid had previously been found for 
P„ = 0fi*2^i22 Ref. m explicitly finds the elastic modes 
for nonzero Pa (recall that a Pq > 25 bars is necessary to 
solidify ^He). A summary of the results and convenient 
notation are provided here. 

For the isotropic case, we define 



a pair that corresponds to first sound and a pair that 
corresponds to fourth sound. With 



f ^ Pi 

Js — , 



(24) 



to first order in /s, first sound frequencies are given by 



fc2 



2c' 



■ w 



(0) 



dp 



1? 



and fourth sound frequencies are given by 



/C2 



„(0) 



dp cl 



Here, the velocities Cq, Ci, and c satisfy 



2 _ dp 



1? 



dP dX 
dp dp 
2 _ dX 



1 /dX dP\ 
p \ dw dw J 



(25) 
(26) 

(27) 
(28) 
(29) 



If a rather than s were held constant, cg would be the 
sound velocity in an ordinary (non-super) liquid, and 
ci would be the sound velocity in an ordinary solid.— 
Ref. shows that for Pa ^ K we have > c? » Cp 

and strain wj^' ^1. It is also convenient to define the 
"fluid- like" and "solid- like" velocities c/l and c/s, which 
satisfy^^ 



-IL 



so that 



dP 
dp 



dX 



I f dX dP 



-IS 



p \dw dw 



(30) 



dP ^dP 
dwik dw 

dX^k r 

dp dp 
dX 4 

d^^'^+r^'- 



(20) 
(21) 
(22) 



In this case the static value of the strain (dependent on 
the applied pressure) is isotropic j^i^i 



w. 



(0) 



w, 



(0) 



3K 



(23) 



1. Longttudinal Elastic Modes 



For k ■ jn ^ k ■ js and A; x {?„ = 0, there are two 
degenerate pairs of solutions to the equations of motion. 



-IS- 



(31) 



For an ordinary solid, the derivatives in (PO]) are taken at 
constant a rather than s. 



2. Transverse Elastic Modes 

For k ■ jn = = k ■ js and k x Vn ^ 0, there are 
two degenerate pairs of elastic modes. They each have a 
frequency satisfying 




(32) 



which is larger than the ordinary (non-super) solid trans- 
verse frequency by the factor yyo/p^. Such an effect, to 
our knowledge, has not been observed. 
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Diffusive Mode 



For the diffusive mode, we keep the dissipative terms in 
the equations of motion ([5|)- (IT^ . so that u[ ^ t>„ ■. With 
ikju[, rewriting ([5])- p^ in terms of the variables 



v'm^ '"U^ P'^ and u'^ gives 

UJS = kiSVni - ifc — ( —S + —p + T^ikjUi I , 



T \ds dp dwji 

bJUi — IVni ~ S H p + — ikiU 



\ ds 



dp dwji 



Wp' = hg[ = ki {pnVn'i + PsVs'i) 



(33) 

(34) 
(35) 



-^Oik I ^ 

OS OS 



, '9P^ dXik \ , 
+ I - — — p 

op op 



dP dXik > ., , 
Oik - I ikjUi 



dw 



(36) 



ujvs'i = hp' = h i^s' + ^p' ' 



ds dp dwji 



tk,u'^ . (37) 



Recall that we have neglected the viscosity as a higher- 
order effect in fc^ as A; ^ 0. We assume that 



'iDok^ 



(38) 



where the diffusion constant Do > is to be determined. 

At first sight this system promises to yield a quintic in 
uj, associated with the longitudinal modes. However, the 
assumption that there is a diffusive mode (whose con- 
sistency we must verify) permits us to reduce this to a 
single linear equation. In some sense a single diffusive 
mode is expected, because we have already obtained four 
pairs of propagating modes. We detail our procedure 
because it both illuminates the physics and clarifies the 
mathematics. 

(1) Method of Solution. Since we take the long 
wavelength limit, we neglect terms that are higher order 
in k. In the present analysis we are merely interested in 
an order of magnitude estimation so we drop subscripts. 
When later solving for the frequency and mode structure 
we use appropriate subscripts. 

When written in terms of powers of k (using ([38])), 
mass and momentum conservation (eqs. psp and (I36p ) 
imply that 



ep' ^kg', 
k^g' r^ka'. 



(39) 
(40) 



Here we use a to distinguish a stress from cr, the en- 
tropy/mass. Combination of ([59]) and pOj) yields k^p' ^ 
a'. Since expanding a' in terms of the other variables 
gives a term proportional to p', for small k the term fc^p' 
is negligible, so cr' as fc 0. The diffusive mode 
therefore is characterized by a negligible stress deviation. 



Physically this means that the fluid-like stress deviation 
nearly cancels the solid-like stress deviation. When a' is 
expanded in terms of the other variables, the condition 
ct' « provides a relationship between s', p' and ku' . 

We now turn to the superfluid equation ((37)) . which 
gives 

fc^l. ~ kp'. (41) 

We now assume that as fc — s- 0, to be verified 

below. When p,' is expanded in terms of the other vari- 
ables, the condition /x' « provides a second relationship 
between s', p' and ku'. In the remaining equations, for 
s' and u', we choose to eliminate p' and u' in favor of s' . 

Neither of the equations for s' or u' (eqs. p3|) - p4)) ') 
involve v'^. Hence, on eliminating p' and u' in favor of 
s', eqs. ([55]) and ([M)) involve s' and v'^, as well as the 
unknown uj. This leaves us with two linear equations for 
two unknowns: the ratio of u'j to s', and w. Once these 
are determined, we use conservation of mass to relate the 
still-unknown v'^ to v'^ and p', both of which having been 
found in terms of s' . We find that at low temperatures 
cjp' can be neglected relative to kv'^, so that ([55)1 gives 
« g' = pnv'^+psv'g. This is not a result of an analysis in 
powers of fc as fc ^ 0, but rather from relations between 
various thermodynamic quantities. 

In what follows, several Maxwell relations from ([T]) are 
used: 

dp^dT dKk^dT_ dX^^_dp_ ^^2) 
ds dp ' ds dwik ' dp dwik 
Further, Ref. 113 gives, for the elastic stress, 

Xik = ~ 2py^ SikWii + pv {wik +Wki) , (43) 

where dX/dw is defined in ((22)) . Since, as in ([23)) . the 

static strain is isotropic (i.e., w[^^ ^ Sik), eq. (j43p im- 
plies that the static elastic stress also is isotropic (i.e., 

'^zk ~ ^ik)- Thus idXik/dp)^.^ and {dXik/ds)^-^ also are 
isotropic, which permits us to define 

dT^^dKk^dX _df^_dKk^dX 
a — — o Oik, „ — ^ — ^ Oik- \^^) 
dwik ds ds dwik dp dp 

Note that eq. gives 

'dw^ = - '^Pv^ SikSji + pv {SijSki + SkjSii) . 

(45) 

(2) Rev^rriting Stress Equation. For the isotropic 
case, using (Iini)-(EI1), & and dH-dUD, eq. dMD 

gives, for negligible total stress. 







fdP 




fdP 


dX\ 


\ds~ 


d-s)'^'^ 


[dp 


" dp) 



Pv — — I ikikiUi — pvik 



dw I 



dd_ 

ds 



hs' + cfj^hp' - {pcfg - pv)ikikiu'i - pvik^u'i, 

(46) 
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where we define 



55 ^ 9A _ aP 
ds ds ds 



(47) 



Since each term of except the last is proportional to 
ki, we have that u'^ is proportional to ki. Thus, kikiu[ = 
k'^u'^, and becomes, on taking the dot product with 
kijk^ and dropping the indices on 







Q(7 I 1 I I'll 

-Q^s +CilP ~ pcigiku . 



(48) 



Further, since u[ ~ ki, substitution of into 
(IMl) gives u^- ^ ki. Then, since (1571) gives ul,- fc^, 
the diffusive mode is purely longitudinal (u^j ^ v'^^ ^ 
u'i ~ ki), and we therefore drop indices for v'^, v[^, and u' 
dotted with k. Moreover, for u'^ ^ ki, (|^5|) gives 



— IK/^KjUi — — — IK Ui. 

owji aw 



(49) 



Note that Yi is independent of T. To evaluate these 
we take Od ~ 25 K,^^ 7714 w 6.7 x IQ-^"^ kg, p sa 2 x 
10^ kg/m^, Pa « 30 bar, and K w 300 bar.^^s Further, 
following evidence that a supersolid phase of ^He can 
only exist at T < 55 mK^ we take T w 50 mK. Then, 
eq. ([56]) yields 



4 4 
m m 

Yi «-4.5 X 10^ ^2 « -3.25 X 10"^ — , 

4 

Fa « -3.25 ^, (57) 



so that Yi ^ Y3 ^ Y2. This inequality applies for any 
T < 55 mK, and therefore applies at any temperature 
relevant to supersolid ^He experiments subject to Pa <C 
K. 

(5) Rewriting s' and u'^ Equations. Substituting 
([5^ . (gHl), and ^ into and into §^ muhiplied by 
—iksYi/Y2, yields 



(3) Rewriting /i' Equation. Since we assume that 
/i' w 0, we neglect the LHS of (l37t : this yields 



w ^s' + ^p' + -^^ikju'i. 
OS op Ow-ji 

Substitution from (l42|) and dH]) gives 

« ^s' + -^p' + ^iku'. 
Op Op op 



(50) 



(51) 



(4) Combining stress and p' equations. Solving 
B8| and (|5T|) for p' and w' gives 



Yi s ' 



-iku' 



Y2S^ 

Y, s' 



(52) 



where we introduce three quantities, each with units of 
velocity to the fourth power: 



dp 



^1 = —^IL '-O^iS' 



9r 9 da 



9p 



p 9s 



p op OS op 



(53) 
(54) 
(55) 



Here we employ (P7| . Ec^. (|5^ holds for any w ~ fc^. 

Appendix |A] uses the results of Ref. to estimate the 
sizes of Yi, Y2, and I3. With the Debye temperature, 
fcs Boltzmann's constant, and 7714 the atomic mass of 
^He, we find 



Yi 



2P^ 



Y. 



Y, 



247r4 keT Pa 
9 9^ rui p ' 
247r4 T^keT K 
9 9^ mi p ' 



(56) 



T \ds sYi dp sFi 
-Y2 '\y2 ds^ Y2 dp dw) ' ■ 



(59) 



We simplify ([58|) - (|59|) by the following argument. We 
take dX/ds to have the same linear T-dependence as 
dP/ds in a harmonic solid, or dX/ds ^ T. Then, since 
s ~ and Y2,Y^ ~ T*, all terms ~ k in dSS]) have 
the same temperature dependence, and the same is true 
for all terms /3 in (|59p . Thus, since Yi ^10® x ¥3 and 
F3 ~ 10 X Y2, in the parentheses of (155)) and (ICTl) the first 
term dominates. Thus (j58p - ([59)) approximately give, on 
rearranging. 



ksv' 



,2KdT 



Y2 ds 



Yo 



(60) 
(61) 



Subtracting (|6T|) from (|60| and dividing by fcs yields 



=ik- 



1 dT 
T~ds 



Y2 dT 



(62) 



which holds for any uj ^ k^. Here we use 
{dX/ds)/{dT/ds) = dX/dT, where p and Wik are implic- 
itly held constant for each derivative. Hence, eqs. ([52)) 
and (|62)) show the ratios of (p', /cu', t;^) to s' to be fre- 
quency independent. 

We now find the frequency of the diffusive mode using 
([SOI) and ([ST)) . Mass conservation from p5l) then relates 
v'g and s', thus yielding all variables in terms of s' . 
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1. Diffusive Mode Frequency 



Cross-multiplication of and (pT|) yields 



.^^KdT _Y2 
T~dl ^ Yi' 



dX 

ds ' 



The frequency of the diffusive mode thus is 

,RdT 



'ik- — , 
T ds ' 



(63) 



(64) 



where 



K + /3sT 



dX/ds 
dT/ds 



1 



Y2 



dx 



(65) 



Here we use Yi 3> Y2. Recall that dX/dT is taken at 
constant p and Wik ■ The frequency thus has a part asso- 
ciated with thermal diffusion (~ k) and a part associated 
with lattice diffusion (^ (3). 

Before finding the full mode structure, it is worth com- 
menting on dMll-dMl)- If l3sT{dX/dT) < k, then we have 
Dd {K/T){dT/ds), as for ordinary thermal diffusion. 
As noted above, however, n is constant in the long wave- 
length limit for the diffusive mode of the supersolid (to 
be verified below), unlike in the case of a fluid or ordi- 
nary solid. Thus, even if the frequency were precisely as 
for normal thermal diffusion, the mode structure (e.g., 
v'n/v'g, etc.) would nonetheless be different than for the 
usual case. 



2. v'„, v's, and u in the Diffusive Mode 
Eq. dSS]) gives 



ijjp_ _ p„v„ 
kps Ps 



(66) 



By ([5^ . the first term on the RHS of is given by 



ojp Lops 



kps kpsS 



Y3 
Yi 



(67) 



Further, using (|55|) and Yi :S>Y2, eq. ((5^ can be written 
as 



Vn' —ik 



R dT s' 
I 

'1 



- 1 



5sr dX 



UJS 

ks 



Y2 R dT 

Yi l3sT dX 
Y2~R~df 



(68) 



Thus the second term on the RHS of ([55]) is given by 



PnV^ upnS 



kpsS 



Yi PsT dX 



Y. 



2 K 



dT 



(69) 



Since experiments^ indicate that p„ > 0.8p, on us- 



ing Yi > ^3, eqs. dSZl) and ^ give -{prj Ps)v'n > 
{iup' /kps). Eq. (|66p therefore becomes, on employing 

una), 



Pn 



-ik 



Pn 



PsS 



1 dT 

T a7 



(70) 



or, equivalently, g' — Pnv'n + Psv'^ ~ 0. Thus the su- 
perfluid velocity is opposite the normal velocity, with a 
weighting given by Pn/ Ps- Since Pnj Ps > 4 we approxi- 
mately have \v'^\ » Note that (|7D|) explicitly relates 
v'g to s', thus completely specifying the eigenmode. 

We now verify the assumption that /i' — >■ for /c — 0. 
We do so by showing the LHS of ([571) to be negligible 
compared to any given term on the RHS (e.g., k'^v'^ <C 
{dp/ds)s'). Counting powers of k, eqs. ((70)) and ((62|) give 
v'^ - -y; - ks'. Thus k'^v'^ - k^s' < {dp/ds)s', which 
shows the consistency of the assumption. 

Furthermore, using eq. ([5^ to write the lattice velocity 
ii' gives 



Comparison to ([68)) yields 



Y2 



Y2UJS^ 

Yi ks ■ 



Yi l3sT dX 
Y2~r' df 



(71) 



(72) 



Since Yi > Y2, unless [PsT{dX/dT)lR] w -Y2/Y1 (an 
unlikely coincidence), we have v'^-^ « — 'u'(Yi/Y2) \u'\. 
Then, by ([70]), for p„ > p^, we have > > 
Since v'^ ^ ii' , mass motion is distinct from lattice 
motion. 



IV. SUMMARY 

We have re-examined the supersolid hydrodynamics of 
Andreev and Lifshitz, including the effects of nonzero 
applied pressure Pa. For Pa 7^ 0, a solid responds with 
both lattice stress Xik and internal pressure P. The de- 
pendence of P and Xik on Pa is found in Ref. [23l . and em- 
ployed here to describe the eigenmodes. We first summa- 
rized the results for the four degenerate pairs of longitu- 
dinal and transverse elastic mode frequencies (including 
fourth sound); because we include Pa and the associated 
strain, the results differ somewhat from those of previous 
work. In addition, again including Pa and the associated 
strain, in the long-wavelength limit we have obtained the 
previously-unstudied diffusive eigenmode. 

The diffusive mode frequency, under certain condi- 
tions, is similar to the frequency of ordinary thermal 
diffusion. However, the mode involves no deviations in 
net stress or net chemical potential, so its properties dif- 
fer from ordinary thermal diffusion. To produce zero 
net stress deviation, the solid-like elasticity component 
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is cancelled by the previously neglected fluid-like compo- 
nent associated with lattice defects. To produce zero net 
chemical potential deviation, the temperature and pres- 
sure deviations must be related. With zero net stress de- 
viation we find that at low temperature there also is zero 
net momentum. With the normal fluid density dominat- 
ing the superfluid density, this means that the superfluid 
velocity is much larger than the normal fluid velocity. 
Because the lattice displacement is coupled to the elas- 
tic strain with a large coefhcient, but the normal fluid 
velocity is coupled to the fluid-like strain (a pressure) 
with a small coefhcient, zero net stress deviation implies 
that the normal fluid velocity is much greater than the 
lattice velocity. This is an unusual phenomenon, since 



in the other modes the lattice velocity and normal fluid 
velocity are nearly equal. 

A previous work studied the lattice diffusion mode 
for a normal solid having distinct velocities associated 
with momentum (e.g., the normal fluid velocity) and lat- 
tice elasticity.— The motivation was to consider that the 
time-delay in the pressurization experiments of Ref. [l^ 
might be due to that mode, under the assumption that 
the sample studied is not supersolid. Similar consider- 
ations can be made for the diffusive mode we have just 
studied, because both modes are diffusive in nature, and 
thus would show a dependence on the sample thickness 
d as (P. A study of this dependence would be of interest, 
to confirm that the effect observed in Ref. [13 is diffusive 
in nature. 
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Appendix A: Relative Sizes of Yi, Y-z and Y3 for 
Small Pa 



Using the results of Refs. 23 and 29 and estimating 
certain thermodynamic derivatives under the condition 
Pa K allows us to estimate the sizes of Yi, Y2, and Y3 . 
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Ref. [23I gives for a normal solid, to lowest order in Pa/K, ((55)) gives, to lowest order in Pa/K, 



,(0) 



Pa 

K 



VPa dK 



dX 

dp "'-^ pK dV 

dP _ y2p2 Q2j^ 

~dp ~ 2pK^ dV^ 

dP 
dw 



-Pa 



V_ OK 

^KdV 



cl 



d\ K 



4 ,, 



(Al) 
(A2) 

(A3) 
(A4) 
(A5) 



Although Ref. '23! evaluates the derivatives of Xik and P 
at constant a rather than s, at T = we have cr « « s, 
so holding either quantity constant should give nearly 
equivalent results for supersolid ^He. All derivatives of 
K here are taken at constant tr, Wiu, and N, so we now 
drop the subscripts. Ref. [H also finds 



Vfl fVd^ dK_ 



(A6) 



Constant Wik constant is equivalent to constant density 
of lattice sites. Because K \s a. measure of the material 
stiffness, one expects K to increase as V decreases, for 
constant N and Wik, i.e., dK/dV < 0. Then by (|X2|) and 
dm we have dX/dp < and dP/dw < 0. 

For the putative supersolid, we approximate da/ds us- 
ing ©, d^, and g^: 



dd_ 

ds 



dX 

'ds 



dP 
ds 



dX dT 



dT 



ds ^ ds ^ dp 



„(o) 



dX 



Recall that, unless otherwise specified, derivatives with 
respect to p, s or Wik are taken with the other two 
held constant. Eq. (IA1[) shows that for Pa ^ we 

have w^^ <C 1. Also,^"^ as noted earlier, p{dT/dp) ~ 
js{dT/ds), where 7 « 10. Thus, 

da dX _i, dT 

d-s'^d-s'^'+^ ^'d-p- 



(A8) 



On neglecting pv, eqs. (|43|) and (|Aip give dX/ds ' 
{dK/ds)w^^^ « {PalK){dK/ds). 

Substitution of dH, dMl, dMl and dH into dSS 



_9A 2 22 

V^Pl \(dK\ 

~ p2^2 \dv) 



2 2^ 
Y2 =S — C,r. + Cq-^ 
p as 



s dX da 



dT 2 



1 

^'d^^W' 
dTK+ \pv 



(A9) 
(AlO) 
(All) 



p dp ds dp dp p 

Note that all terms ~ da/ds are higher order in Pa/K 
and therefore are neglected. Approximating K to be lin- 
ear in V and neglecting pvi eqs. dA9P - dAlip give 



2^ 

r.2 ' 



dTPa 
dp p 



dT K 
'dp p 



(A12) 



For Pa ^ K we have Y3 > 1^2- 

To approximate the relative magnitudes of ¥3 and Y2 
to Yi, we now find an explicit form for s{dT/dp)s- 

At low temperatures phonon gas statistical mechanics 



gives 



27r24T3 



(A13) 



where u is an average sound velocity and fc^ is the Boltz- 
mann constant. Further,— 



dT ^Tdu ^ lOT 
dp u dp 3/9 

Combining dA13p and (IA14P gives 



dT _ 47r2fc4 



dp 9ph^u^ 



(A14) 



(A15) 



In terms of the Debye temperature 9d ~ 
(67r^ni/)^/"^(te/fcB), where is the number density of 
vibrations (essentially one per lattice site). 



dT 247r4fcBr4 247r'* keT 



dp 9{p/n^ 



9 el 



(A16) 



Here 777.4 is the atomic mass of ^He, and we have taken 
777477^ w p. Eq. dAiep substituted into dA12p gives d56)) . 



